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. Abstract. We study hadron production in jets by applying quark recombination to jet shower 

r-{ ' partons. With the jet showers obtained from PYTHIA and augmented by additional non- 

perturbative effects, we compute hadron spectra in e + + e~ collisions at = 200 GeV. 
Including contributions from resonance decays, we find that the resulting transverse momentum 
spectra for pions, kaons, and protons reproduce reasonably those from the string fragmentation 
as implemented in PYTHIA. 
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. 1. Introduction 

Hadron production from energetic jets in high-energy nuclear collisions is often parameterized 
CN ■ through fragmentation functions using the universality of the process as given by QCD 

factorization theorems [Tj. On a more microscopic level, one can achieve a good description 
of jet hadronization through a perturbative evolution of the parton spectrum inside the jet 
using DGLAP splitting kernels from the largest momentum scale to some lower cutoff Qo, 
| followed by a non-perturbative hadronization model like the Lund string model applied to the 

parton shower. The event generator PYTHIA j2] has successfully implemented such a strategy 
to describe hadron production in e + + e~ collisions at high energies and in numerous other 
processes. 

On the other hand, it has been suggested that hadron production can often be described 
through a process of quark recombination or coalescence [H 12 [5j EJ [7] . It is an intriguing idea to 
combine the concepts of quark recombination and parton showers since it would easily generalize 
to the hadronization of jets in dense environments as in relativistic heavy ion collisions. Hwa 
and Yang have pioneered the application of quark recombination to jet showers in Ref. [8]. 
However, the parton showers in their work were not obtained from first principles but fitted to 
measured hadron spectra. They also ignored the smallness of the number of partons in each jet 
by applying recombination to averaged spectra instead of event-by-event showers. Fluctuations 
from shower to shower are expected to lead to significant changes to the hadron spectrum at 
large hadron momentum fraction z. 

Using e + + e~ collisions as an example, we show that vacuum fragmentation functions can 
be reproduced in an event-by-event approach that uses a perturbative shower evolution in the 
vacuum down to a scale Qo, a minimal modeling of non-perturbative effects on partons below 
Qo, and finally recombination of the quarks and antiquarks into hadrons. 



2. Shower Partons from Electron-Positron Annihilation in PYTHIA 
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Figure 1. (Color online) Transverse momentum spectra of shower partons before (left panel) 
and after (right panel) including gluon decays into quark-antiquark pairs. 



Using PYTHIA 6.3 [2], we generate an ensemble of parton showers from e + + e~ collisions 
at center-of-mass energy y/s = 200 GeV. The light quark, strange quark and gluon spectra as 
functions of momentum at a cutoff scale Qq = 1 GeV are shown in the left panel of Fig. [TJ 
Since gluons in the parton showers cannot directly recombine into hadrons, we force them to 
decay into quark-antiquark pairs after the perturbative evolution is stopped. This is a non- 
perturbative process, and we use constituent quark masses of m u d = 0.3 GeV for light quarks 
and m s = 0.45 GeV for strange quarks with the remaining gluon virtuality randomly chosen in 
the interval between 2m u ^ and Qq. The transverse momentum spectra of the decay light and 
strange quarks together with the total light and strange quark spectra are shown in the right 
panel of Fig. [TJ 



3. Quark Recombination 

Following the approach to quark recombination in Ref. [5] the probability for quarks and 
antiquarks to form a hadron is taken to be proportional to the quark Wigner function of the 
hadron. Assuming the quark wave functions of the hadrons to be those of harmonic oscillators, 
the hadron Wigner function is then Gaussian in both the coordinate and the momentum space. 
Here we neglect spatial information, assuming shower partons with similar momenta to be close 
in space. Hence we only consider the momentum-space Wigner function. 

For the production of a meson from a quark and an antiquark, the probabilities are 

/ k 2 \ / 16 k 2 \ / k 2 \ 

/s(k) = 8fif M exp ( --^2 j and / p (k) = g M ( -8 + j exp ( -^2 j , (1) 

respectively, for s-wave and p-wave mesons [9]. Here k is the relative momentum between 
the quark and antiquark in the meson rest frame, A is the width of the quark and antiquark 
relative momentum distribution in the meson, and is the statistical factor for the colored, 
spin-1/2 quark and antiquark to form a colorless meson. For the mesons tt, p, ai, U), /2(1270), 
/l(1285), K*, 2^(1270), and XJ(1410) considered in our study, we have g n = g K = 1/36, 
g P = 9u = g ai = 9f 1 = 9K* = = 9K* = 1/12, and g h = 5/36. Because of the color charge 



of a gluon, its decay quark and antiquark are excluded from recombining into a color singlet 
meson. 

For s-wave baryons, the production probability from three quarks is [10] 

/ k 2 Sk 2 \ 

/ fl (k! , k 2 ) = Ug B exp ( - - |^) , (2) 

where ki = (p^ — p 2 )/2 and k 2 = (p^ + p 2 — 2p(j) /3 are the relative momenta of the recombining 
quarks of momenta p^, p 2 , and P3 in the rest frame of the baryon. The statistical factors for 
the baryons N, A?i(1440), A, and Ai(1600) considered in our study are g^ = g^ = 1/108 and 
9A = 9Ai = 1/54. 

4. Pion, Kaon and Nucleon Transverse Momentum Spectra 
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Figure 2. (Color online) Transverse momentum spectra of pions (left panel), kaons (middle 
panel), and nucleons (right panel). 

The momentum width A is related to the oscillator constant in the quark wave function. 
Here we take them as parameters. Using A^ = A p = A ai = A w = Af 2 = Af 1 = 0.43 GeV, 
Ak = Ak* = Ak ± = Ak* = 0.58 GeV, and Ajv = Ajv x = Aai = 0.62 GeV, we have calculated 
their respective transverse momentum spectra from the recombination of quarks and antiquarks 
on an event-by-event basis. The spectra of directly produced pions, kaons and nucleons are shown 
by star symbols in the left, middle, and right panels of Fig. O respectively. They underestimate 
those obtained from PYTHIA via string fragmentation shown by solid lines. After including 
contributions from the decay of resonances via 

p — > 7T7T, dl — > pTT — > TTTTTT, UJ —> 7T7T7r , / 2 — > 7T7T , f\ — > 4lT , 

K* -> Kit, K x -> Kp -> Ktttt, K{ -> K*tt -> Ktttt, 
N t -> A%(65%) or Np{35%) -»• A%vr, 
A -)■ Ntt, Ai -)• Ntt, 

we significantly improve the comparison with the PYTHIA results as shown by the dashed lines 
in Fig. [2 

We note that not all quarks and antiquarks are used in the recombination model to form 
hadrons. As shown in Fig. [3l the recombination probability is close to one for low momentum 
quarks and decreases with increasing momentum. The small number of non-recombined quarks 
thus consists mostly of leading quarks that are far removed in momentum space from the rest 
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Figure 3. (Color online) Probabilities of light and strange quarks to recombine into hadrons 
as functions of their transverse momenta. 



of the parton shower. Their contribution to hadron production can be included via independent 
fragmentations at the low momentum scale Qq. Because of their small numbers, the final hadron 
spectra are, however, not much affected. 

5. Conclusions 

Our study indicates that hadron production from jets can be described by the recombination 
of shower partons inside the jets. Our approach provides a promising framework to study the 
medium modification of fragmentation functions in the presence of a quark gluon plasma (QGP) 
by allowing shower partons to recombine with the thermal partons in the QGP. 
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